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I.  INTRODUCTION 
1.1   Background 

Although  extensive  ocean  ambient  noise  studies  were  conducted 
during  World  War  II  with  continuing  effort  since  that  date,  few  of 
either  theoretical  or  experimental  treatments  hare  dealt  with  the 
directional  properties  of  such  noise.  There  are  many  possible  sour- 
ces of  ambient  noise  in  the  ocean  any  of  which  could  produce  non- 
uniform distributions.  R.  A.  Frosch^  '  in  examining  likely  contrib- 
utors to  the  total  acoustic  energy  content  of  the  oceans  has  con- 
sidered, for  example,  biological  sources,  seismic  noise,  shipping, 
lightning  and  sea  surface  waves.  His  conclusion  that  sea  surface 
waves  are  a  major  contributor  is  in  agreement  with  World  War  II 
data.  For  example,  the  "Knudsen"'2'  curves,  a  compilation  of  omni- 
directional measurements  relating  ambient  noise  spectrum  levels  to 
frequency  as  a  function  of  sea  state,  exhibit  a  direct  sea  state/ 
ambient  noise  interdependence  and  by  their  very  nature  permit  the 
inference  of  noise  anisotropy  in  the  ocean. 

Theoretical  treatments  considering  the  ocean  surface  as  a 
possible  noise  source  have  been  published  by  Eekart^*',  Roberson^  ', 
Urick^"'  and  Miller'  '.  Eckart  has  presented  the  physical  picture  of 
the  ocean  surfaoe  resolved  into  sinusoidal  travelling  waves  of  the 
particular  phase  velocity  which  would  be  associated  with  free  flex- 
ural  waves  of  an  elastic  plate  excited  by  random  local  shocks  of 
short  duration.  The  above  model  would  result  in  the  generation  of 
sound  only  under  the  unlikely  condition  that  the  surface  phase 
velocity  exceeded  the  velocity  of  the  bounded  medium.  Roberson 
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without  attention  to  a  phyaioal  model  of  noise  production  considered 
both  an  infinite  volume  distribution  of  simple  sources  and  the  com- 
bined effect  of  a  volume  distribution  of  simple  sources  bounded  by  a 
surface  distribution  of  dipole  sources*   Assuming  an  attenuation  law 
proportional  to  frequency  squared,  he  showed  a  -6  decibel/octave  de- 
pendence of  noise  upon  frequency  for  the  unbounded  case  and  a  non- 
uniform dependence  upon  frequency  for  the  mixed  case  depending  upon 
the  ratio  of  the  contributions  from  the  two  sources  and  the  fre- 
quency. In  this  latter  case,~'no  combination  of  surfaoe  and  volume 
sources  produced  a  constant  slope  in  agreement  with  the  "Knudeen" 
curves  except  over  a  limited  frequency  range. 

Urick  specifically  considered  the  subject  of  directional 
distribution  of  ambient  noise  and  concluded  that  significant  aniso- 
tropy  might  exist  depending  upon  the  unknown  manner  in  which  each 
element  of  surface  radiates. 

Miller  considered  the  directional  distribution  resulting 
from  the  combination  of  straight  path  propagation  from  surfaoe  rad- 
iators close  to  the  receiver  and  convergent  zone  propagation  from 
distant  surface  radiators.  The  acoustic  radiation  from  the  surface 
varied  as  cosine  8  ,  where  8   is  the  angle  between  the  direction  of 
radiation  and  a  normal  to  the  radiating  surfaoe.  No  bottom  reflec- 
tion or  volume  sources  were  assumed.  He  concluded,  tentatively, 
that  noise  in  the  above  case  is  essentially  isotropic  over  a  hemis- 
phere with  a  slight  increase  within  ♦  9  degrees  of  the  horizontal 
and  that  the  area  immediately  above  and  close  to  the  reoeiver 


contributes  about  the  sane  total  noise  as  the  contribution  from  the 
first  convergence  zone. 

If  the  ocean  surface  is  a  major  contributor  to  ambient  noise, 
it  has  been  reasoned  that  ambient  level  should  decrease  with  depth. 
The  change  in  ambient  noise  levels  with  deptrn^)  has  been  investi- 
gated experimentally.  No  such  change  has  been  detected  to  depths  of 
500  feet  in  tests  at  25  kilocycles  per  second.  Urick'5)  has  shown, 
however,  that  the  change  to  be  expected  between  5  a**d  500  feet  depth, 
even  at  a  higher  frequency  of~50  kilocycles  per  seoond  where  attenu- 
ation is  greater,  is  only  5  decibels  theoretically  and  thus  a  smaller 
change  would  have  been  obscured  by  measurement  fluctuations  at  the 
lower  frequencies.  Teer^"/  reported  the  results  of  tests  at  19 
kilocycles  per  second  with  a  small  line  array  oriented  horizontally 
and  vertically  in  deep  water.  The  results,  which  seemed  to  indicate 
that  ambient  noise  is  isotropic,  are  not  considered  conclusive  be- 
oause  of  limitations  in  equipment  and  the  environmental  conditions 
experienced  during  the  tests. 

In  summary,  no  experimental  data  exist  to  describe  the  di- 
rectional distribution  of  ocean  ambient  noise.  Theoretical  treat- 
ments have  indicated  either  isotropic  or  anisotropic  distributions 
may  exist  depending  upon  initial  assumptions  as  to  the  nature  of  the 
surface  radiator  and  the  uncertain  effect  of  ocean  sound  channels. 

1 .2  Research  Objectives 

It  is  the  purpose  of  this  study  to  investigate  anisotropy 
of  deep  ocean  ambient  noise  in  the  frequency  range  of  750  to  1500 


cycles  per  second.     Further  specific   objectives  will  bet 

a.  To  isolate  vertical   and  azimuthal   dependencies. 

b.  To  establish  the  identity  of  the  sources  contributing 
to  non-uniform  distributions. 

c.  To  infer  the  radiation  patterns  of  such  sourees. 


II,   AMBIENT  NOISE  EQUIPMENT 
2.1   Discussion 

There  are  two  major  reasons  which  hare  prevented  more  wide- 
spread Investigation  of  the  directional  properties  of  ambient  noise. 
The  first  concerns  the  limitation  of  a  single  directional  hydrophone 
which,  to  monitor  several  directions  in  space  successively,  requires 
the  assumption  of  a  stationary  process.  The  above  difficulty  may  be 
overcome  by  the  employment  of  more  than  one  directional  hydrophone 
or  by  forming  several  directional  beams  from  an  array  of  fixed  non- 
directional  hydrophone  elements.  Either  solution  generates  a  second 
problem,  namely  the  difficulty  of  ensuring  uniform  gain  in  the  sev- 
eral beams.  Both  limitations  have  been  surmounted  by  a  system  oon- 
ceived  by  Anderson^"'  and  implemented  by  the  staff  of  the  Marine 
Physieal  Laboratory  of  the  Serippc  Institution  of  Oceanography, 
University  of  California,  La  Jolla.  A  52-element  array  has  been 
constructed  in  the  shape  of  a  great  stellated  loosohedron  as  shown 
in  Fig.  1.  A  summary  of  the  array  design  parameters  is  eontained  in 
paragraph  2.2  following.  The  52  hydrophone  outputs  are  combined  to 
form  ^2  beams  with  axes  coinciding  with  lines  passing  through  the 
array  center  and  the  individual  hydrophones.  All  52  hydrophones  are 
used  in  forming  each  beam. 

The  problem  of  uniform  sensitivity  in  each  beam  is  minimized 
by  employing  clipped  processing.  At  the  output  of  eaoh  hydrophone 
channel  prior  to  beam  formation,  all  amplitude  information  is  dis- 
carded by  infinite  clipping  and  phase  information  only  retained. 
The  clipped  "on-off"  form  of  the  channel  outputs  permits  digital 


FIGURE 


THE    AMBIENT    NOISE    ARRAY 


shift  register  techniques  to  be  employed  in  the  forming  of  the  52 
beams  with  an  attendant  faring  in  space  and  equipment  complexity 
over  conventional  analog  beam  forming  techniques.  The  system  is 
thus  independent  of  individual  ohannel  gain  differences  from  each 
hydrophone  through  the  beam  forming  process*  The  problem  of  ampli- 
tude differences  appears  only  in  compensating  for  detector  gain  dif- 
ferences in  the  process  of  rectifying  each  beam  output.  The  dipping 
and  beam  forming  processes  utilised  in  this  system  have  been  re- 
ported by  Anderson' 10'  and  Rudnick'   '.  The  restoration  of  the  ana- 
log form  of  a  signal  from  clipped  hydrophone  outputs  is  considered 
briefly  in  paragraph  2.% 

2.2  Array  Configuration  Development 

In  arriving  at  the  present  array  geometry  of  a  great  stel- 
lated ieosohedron.  Anderson'"'  considered  three  factors! 

a.  The  number  of  elements* 

b.  Element  configuration* 

c.  Element  spacing. 

Since  presumably  small  deviations  from  an  isotropic  field  were  to  be 
measured,  a  symmetrical  distribution  was  suggested  permitting  nearly 
identical  directivity  patterns  of  the  several  beams  to  be  formed  by 
electrical  steering.  A  distribution  of  elements  in  depth  rather 
than  a  spherical  shell  would  permit  greater  interelement  spaeings 
for  a  given  array  diameter  and  thus  increased  array  efficiency 
through  minimising  interelement  correlations.  The  maaber  of  elements 
was  limited  by  considerations  of  electronic  complexity.  Thirty-two 
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appeared  a  practical  number  and  (a)  and  (b)  above  were  satisfied  by 
the  great  stellated  icosohedron.  Optimum  design  of  intorolement 
spacing  was  determined  by  calculations  of  array  efficiency  defined 
as  the  ratio  of  output  to  input  signal  to  noise  ratios.   Signal  in 
this  case  was  the  deviation  from  isotropy  which  it  was  desired  to 
measure  in  a  given  direction.  An  interelement  spacing  of  0.6  wave 
lengths  at  the  1  kilocycle  per  second  nominal  center  frequency  was 
chosen. 

2.3  General  System  Description 

An  overall  system  block  diagram  is  shown  in  Fig.  2.  For 
simplicity  only  1  of  52  identical  receiving  channels  is  represented 
completely.  Reading  from  the  left,  an  individual  hydrophone  output 
is  first  amplified  and  then  employed  to  frequency  modulate  a  tele- 
metering transmitter.  The  52  transmitter  carriers  are  distributed 
from  10  through  16.2  megacycles  per  second  in  200  kilocycles  per 
second  increments.  An  additional  1  megacycle  per   second  amplitude 
modulated  transmitter  relays  hydrophone  array  asdmuthal  orientation 
in  accordance  with  the  directional  indication  of  a  magnetic  compass. 
The  52  hydrophones,  preamplifiers,  telemetering  transmitters  and 
compass  system  comprise  the  underwater  portion  of  the  system. 

The  transmitter  outputs  are  coupled  to  the  shipborne  re- 
ceiving portion  by  coaxial  cable.  A  1500  foot  length  was  available. 
The  5/8  inch  diameter,  doubly  armored  cable  serves  also  as  the 
supporting  member  of  the  700  pound  receiving  array.  Cable  attenua- 
tions of  1.5  and  1.8  decibels  per  100  feet  at  10  and  16  megacycles 
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FIGURE  2.   AMBIENT  NOISE   EQUIPMENT    BLOCK    DIAGRAM 
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per  second  respectively  necessitated  the  insertion  of  a  broadband  10 
to  1^,2  megacycle  per  second  preamplifier  prior  to  demodulation  of 
the  signal  by  the  "$2  matching  frequency  modulated  receivers.  At  thia 
stage  it  iR  possible  to  view  the  amplitude  and  spacing  of  the  J>2  in- 
dividual carriers  as  v/ell  as  the  frequency  deviation  imposed  by  each 
signal  by  means  of  a  panoramic  receiver  and  cathode  ray  tube  indica- 
tor. 

Following  demodulation  in  the  FM  receivers,  the  audio  signal 
is  further  amplified  and  then  infinitely  clipped.  All  amplitude  in- 
formation is  discarded  and  only  polarity  information  retained.  The 
clipped  signal  then  enters  a  digital  shift  register  ten  bit  delay 
line  which  provides  suitable  signal  delay  increments  such  that  sub- 
sequent combination  of  the  J>2   signals  in  a  mixing  matrix  yields  52 
beams  symmetrically  spaced  over  4  zd2.  The  J>2   beam  outputs  are  then 
individually  detected  and  integrated  by  a  simple  low  pass  RC  filter 
of  one  second  time  constant.  Following  detection  and  integration, 
the  p2  beam  outputs  are  sequentially  sampled  at  a  50  cycle  per  second 
rate  by  a  mercury  ^et  scanning  switch  for  presentation  on  a  cathode 
ray  oscilloscope.  Sixteen  millimeter  scope  photographs  preserve  the 
data  for  analysis. 

All  receiving  equipment  is  housed  in  a  portable  Army  Signal 
Corps  communications  hut.  The  entire  7,000  pound  installation  may 
be  moved  quickly  by  crane  from  its  laboratory  position  to  the  fantail 
of  PA0LINA  T,  a  Scripps  Institution  of  Oceanography  research  vessel. 
The  electronic  hut,  array,  associated  winches  and  boom  and  gasoline 
powered  15  KW  generator  satisfying  the  115  volt,  50  ampere  demand 
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are   shown  in  Fige.   5  and  4. 

2.4     Hydrophone  and  Transmitter  Assembly 

The  hydrophones  are  encapsulated  barium  titanate  cylinders 
with  a  flat  response  averaging  90  decibels  below  1  volt//^    bar 
throughout  the  range  of  interest,   750-1500  cycles  per  second*     The 
cylindrical   sensitive  element  is  1-1/2  inches  high,   1-7/16  inches  in 
diameter  and  has  a  wall  thickness  of  1/8  inch.     Individual  hydro- 
phone directivity  is  discussed  in  paragraph  5*5.     Although  the  par- 
ticular units  employed  were  constructed  at  the  Marine  Physical   Lab- 
oratory,  an  identical  unit,   the  REMACO  Model   R-1J0  is  commercially 
available.     They  are  rated  at  a  working  pressure  in  excess  of  1600 
pounds  per  square  inch.     The  hydrophone  fits  into  the  end  of  a 
16-1/2  inch  long,  5/52  inch  thick  stainless  steel  pipe  which  forms 
the  watertight  housing  for  preamplifier  and  transmitter.     All  hydro- 
phone units  are  series  connected,   50  VDC  required  by  each  unit  or 
1600  VDC  total.     A  relay  feature  is  included  to  remove  a  particular 
hydrophone  unit  should  a  vacuum  tube  filament  open.     Were  this  fea- 
ture not  included  the  entire  array  would  beoome  disabled  by  the 
failure  of  a  single  tube.     On  the  other  hand,  a  short  circuit  any- 
where in  the  series  string  disables  all  units  beyond  the  point  of 
shorting.     A  partially  disassembled  unit  is  shown  in  Pig.   5» 

2«;i     Clipped  Processing 

Because  digital  multibeam  steering   (DIMUS)   is  a  relatively 
new  technique,   the  following  simplified  illustration  of  the  restora- 
tion of  the  analog  form  of  a  signal  by  the  addition  of  several 
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samples  of  the  signel  which  have  been  clipped  is  included.  Referring 
to  Fig,  6a,  the  time  varying  analog  audio  output  of  an  individual 
hydrophone  channel  is  shown.   Some  band  limiting  has  occurred  in  the 
hydrophone  preamplifier  bandpass  filters.  Further  high-pass  filter- 
ing by  the  simple  network  of  Fig.  6b  removes  any  remaining  fluctua- 
tion below  the  frequency  of  interest  resulting  in  the  waveform  of 
Fig.  6c.  The  waveform  of  6c  is  transformed  into  one  of  two  states 
as  shown  in  Fig.  6d  by  the  dipper  amplifiers.  Figure  7a  repeats 
Fig.  6d.  In  addition  Figs.  7b, ^c,  d  and  e  represent  the  outputs  of 
four  other  channels.  The  example  has  been  so  chosen  as  to  illustrate 
the  array  gain  of  a  five  hydrophone  system  in  detecting  a  directional 
signal  mixed  with  omnidirectional  noise.  An  input  signal -to-noise 
ratio  at  the  individual  hydrophone  of  unity  has  been  assumed.  The 
direct  addition  of  the  five  output  voltages  produces  a  randomly 
fluctuating  voltage  as  shown  in  Fig.  7f •  The  absence  of  a  pronounced 
peak  is  to  be  expected  since  the  noise  at  each  hydrophone  due  to  the 
isotropic  background  is  uncorrelated  while  correlation  between  volt- 
ages generated  by  the  directional  noise  will  be  visible  only  when 
the  proper  electrical  delays  corresponding  to  the  geometrical  delays 
between  hydrophones  are  introduced  by  the  shift  registers.  The 
voltages  of  Figs.  7a  through  c  are  sampled  at  the  shift  register  in- 
puts at  a  5  kilocycle  per  second  rate.  This  timing  satisfies  the 
Nyquist  criterion  of  a  minimum  sampling  rate  equal  to  twice  the 
highest  frequency  of  the  spectrum  of  interest.  Figures  7g  through 
7k  repeat  Figs.  7a  through  7e  with  the  proper  delay  introduced  by 
the  shift  registers  to  correct  for  hydrophone  geometry.  The 
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addition  of  Figs.  7g  through  7k  shown  in  Fig,  71  again  produces  a 
randomly  fluctuating  voltage  with  the  exception  that  now  a  peak  ap- 
pears due  to  the  directional  noise.  The  reproduction  of  the  original 
signal  improves  with  the  number  of  channels  employed.   In  aural  test 
with  the  52  element  system  it  is  not  possible  to  distinguish  between 
an  undipped  audio  signal  input  as  monitored  at  the  output  of  a  sin- 
gle hydrophone  and  the  reconstituted  input  as  heard  at  the  output  of 
a  beam  formed  from  the  clipped  signals  of  the  32  elements. 

2.6  Array  Beam  Patterns 

The  summed  voltage  output  of  all  hydrophones  may  be  written: 

%ut--£z'^L   fcos(^-^.)cos(^-r</)-cos(4-^  +  ^)cos(4-^+^)' 

is  I  c  L  J 

where        cj     ■  angular  frequency,    radians  per  second 
RnOl,<pi    m  hydrophone  coordinates 

&s><fis   *  electrically  steered  direction 
&  ,<fi     m  angular  deviations  from  steered  direction 
Due  to  hydrophone  symmetry,   a  particular  beam  will   exhibit  conical 
symmetry  about  the  beam  axis.     The  icosohedron  has  two  axes  of  sym- 
metry and,   accordingly,   two  different  beam  patterns.     In  a   spherioal 
coordinate   system  with  the  array  oenter  as  origin,   every  element  can 
be  pictured  as  positioned  on  one  of  two  concentric   shells.  There  are 
20  elements  located  on  the  outer  shell  and  12  elements  on  the  inner 
shell.     The  two  axes  of  symmetry  are  formed  by  lines  through  the 
array  center  and  elements  on  the  outer  and  inner  shells.     Since  each 
array  element  has  a  diametrically  located  twin,  the  expression  for 
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power  output  with   8S   -    <fis  =  0     and     <f>D    =  0     becomes: 

G{&D)=J     |2^cos(f^-/    cos^.  cos^-cos(<^-<?)    )}  B{w)dou 

where  5* (6;)  i3  a  weighting  function,  the  relative  frequency  response 
of  the  system  and  ^ ,  cu^    are  the  system  bandpass  limits.   Computed 
patterns  normalized  to  unit  on-axis  response  are  shown  in  Fig.  8. 
The  continuous  curve  represents  the  solution  to  the  preoeding  equa- 
tion for  an  outer  beam  (a  beam  whose  axis  passes  through  an  outer 
element  in  the  array). 


The  use  of  the  preceding  equation  as  written  tacitly  assumes 

that  the  exact  time  delay  specified  by  the  expression  -cos$2>  COS<# 

it  available.  In  practice  the  exact  delays  are  not  available  owing 


to  the  limitations  of  the  digital  shift  register  which  at  its  tapped 
outputs  provides  delays  equal  to  integer  multiples  of  the  clock  gen- 
erator period.  The  dashed  curve  in  Fig.  8  is  the  outer  beam  re- 
sponse computed  on  the  basis  of  the  actual  shift  register  delays. 
The  dotted  curve  Is  the  corresponding  response  for  an  inner  beam. 
The  combined  effect  of  beam  differences  due  to  element  configura- 
tions and  steering  delays  results  in  an  inner  to  outer  on-axis  beam 
response  ratio  of  O.855,  a  factor  which  must  be  considered  in  the 
data  reduction  process. 

2.7  Reduction  of  External  Noise  Interference 

The  preliminary  phases  of  at  sea  measurements  confirmed  the 
existanee  of  two  sources  of  interfering  noise  which  had  been  antici- 
pated. They  were: 
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a.  Ship  machinery  noise  transmitted  through  the  hull  struc- 
ture into  the  ocean. 

b.  Self-noise  generated  by  relative  movement  of  the  array 
and  ocean  water. 

Another  troublesome  source,  external  shipping  interference,  eould 
only  be  controlled  by  Judicious  choice  of  operating  area  free  from 
shipping  lanes  or  by  testing  during  night  time  hours  when  shipping 
traffic  was  at  a  minimum.   A  fourth  source  generated  by  the  rolling 
motion  of  the  testing  ship  is  discussed  in  paragraph  %1« 

Own  ship  machinery  noise  was  eliminated  by  two  means.  First, 
all  normal  ship  machinery  was  secured  during  the  testing  period. 
Secondly,  interference  from  the  speoial  15  KW  gasoline  powered  gen- 
erator serving  as  the  power  supply  for  the  ambient  noise  equipment 
was  reduced  to  a  negligible  value  by  isolating  the  generator  from 
the  ship  deck  by  automobile  tire  casings.  The  effectiveness  of  the 
isolation  method  was  demonstrated  by  comparative  evaluation  with  a 
special  5  KW  sound  isolated  unit  developed  at  the  U.  S.  Navy  Elec- 
tronics Laboratory  for  similar  noise  teste.  This  latter  source  was 
completely  suspended  from  each  of  its  four  corners  in  a  pipe  support 
frame  by  airplane  shock  cord.  Although  its  electrical  capacity  was 
not  adequate  for  full  time  employment,  alternate  operation  with 
first  the  limited  capacity  *  KW  unit  and  then  the  larger  15  KW 
source  revealed  no  differences  in  ambient  noise  levels  as  measured 
at  the  array  suspended  in  the  ocean. 

Elimination  of  the  second  interfering  source,  self-noise 
generated  by  relative  array/water  motion,  was  accomplished  by  two 
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means:   (1)  the  suspension  of  the  array  from  a  system  of  buoys  and 
floats  in  order  to  isolate  the  array  from  ocean  surface  and  ship  mo- 
tion and  (2)  the  employment  of  a  parachute  sea  anchor  to  eliminate 
Bhip  drift  caused  by  surface  winds  and  attendant  array/water  rela- 
tive motion.  The  two  techniques  are  shown  schematically  In  Fig.  9« 
The  suspension  was  so  designed  that  the  spar  buoy  would  float  in  the 
position  shown  independent  of  array  depth.  The  floats  compensated 
for  cable  weight  as  the  array  was  positioned  at  100,  500  or  1000  foot 
depths.  The  spar  buoy  acts  effectively  as  a  spring  with  constant, 
k  ,  equal  to  the  buoyancy  force  per  length  submerged.  For  the  4  Inch 
diameter  buoy  in  question  k  »  5  lbs/ sec  which  for  an  array  dis- 
placement of  500  pounds  results  in  a  natural  period  of  oscillation  of 

T  -  7.ir^/^r-   ~  60  seconds 
V  k 

a  value  10  times  greater  than  the  longest  period  of  the  ocean  waves 
encountered  during  the  experiment.  Particular  attention  was  devoted 
to  buoy  swivels  and  connecting  links  in  order  to  prevent  noise  gen- 
eration by  any  metal -to-metal  contact. 

The  sea  anchor  employed  was  a  surplus  government  parachute 
fastened  in  the  manner  shown  by  manila  hawser  and  equipped  with  a 
manila  tripping  line  for  collapsing  purposes.  Even  in  very  light 
wind  of  7  knots  it  was  effective  in  maintaining  ship's  head  into  the 
seas,  reducing  ship  roll  and  completely  eliminating  drift.  Under 
certain  not  infrequent  conditions  stratified  ocean  currents  were  ex- 
perienced where  surface  currents  and  wind  were  in  a  direction  differ- 
ent from  the  ocean  current  at  the  array  depth.  The  net  effect  was  "the 
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tovlng  of  the  array  by  the  teat  ship  and  parachute  and  invariably 
resulted  in  test  cancellation. 


III.   SYSTEM  MEASUREMENT  ACCURACY 
%1   Discussion 

The  ideal  manner  in  vhieh  to  specify  system  measurement  ac- 
curacy would  be  in  terms  of  beam  response  to  a  calibrated  source  cap- 
able of  producing  an  output  spectrum  covering  the  ambient  noise 
equipment  bandwidth.  Up  to  the  present  the  above  procedure  has  not 
been  possible  for  a  variety  of  reasons.  The  chief  limitation  arises 
from  the  low  frequency  of  operation  of  the  system.  Neither  the  nec- 
essary low  frequency  broadband  source  nor  a  quiet,  deep  water  area 
with  adequate  facilities  for  maintaining  array/source  geometry  has 
been  available  for  test.  Although  it  has  not  been  possible  to  cali- 
brate the  system  as  a  single  unit  by  one  series  of  measurements, 
calibration  of  all  component  parts  of  the  system,  individually,  has 
been  possible  and  is  considered  an  adequate  substitution  for  the 
single  measurement  procedure.   In  paragraphs  %2  through  3«7  follow- 
ing, the  possible  error  contributions  of  all  component  parts  of  the 
system  are  examined.  Deteotor  calibration  has  been  established  as 
the  major  sOuree  of  error  in  the  system  and  will  accordingly  be  con- 
sidered first.   It  is  the  objective  of  the  remaining  paragraphs  to 
show  that  all  contributions  from  other  components  of  the  system  are 
negligible  in  contrast  to  deteotor  calibration  error. 

3»2  Detector  Calibration 

The  use  of  clipped  processing  considerably  complicated  de- 
tector calibration  prooedures.   Since  the  output  amplitude  of  a 
given  beam  is  a  function  of  the  relative  phase  correlation  between 
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individual  hydrophones  and  is  independent  of  individual  hydrophone 
output  amplitude,  the  input/ out put  characteristic  of  a  given  beam 
forming/detector  combination  could  not  be  established  by  injecting  a 
signal  of  controlled  amplitude  into  each  dipper  amplifier  and  mea- 
suring detector  output  voltage. 

The  problem  is  further  complicated  by  the  nature  of  the  shift 
register  matrix  output.  A  representative  waveform  in  shown  in  Pig, 
10a  illustrating  the  summed  voltage  output  from  52  shift  registers 
with  random  noise  input.  The  desired  information  is  contained  in 
the  fluctuating  AC.  The  undesired  5  kilocycle  per  second  clock 
pulses,  arising  from  the  particular  form  of  shift  register  circuitry 
employed,  eliminated  the  possibility  of  employing  full  wave  rectifi- 
cation and  necessitated  use  of  less  efficient  half-wave  rectifiers. 
A  quantitative  analysis  of  the  system  output  including  the  effects 
of  clipping  is  contained  in  Appendix  I. 

Although  52  random  noise  generators  whose  relative  correla- 
tions could  be  varied  at  will  would  permit  detector  calibration,  no 
circuitry  to  fulfill  this  requirement  was  obvious.  Instead  a  per- 
fectly correlated  AC  signal  was  employed.  A  2.5  kilocycle  per  second 
test  voltage  of  two  phases,  180  degrees  apart,  was  derived  by  fre- 
quency divider  from  and  synchronised  with  the  5  kilocycle  per  second 
clock  generator.  Application  of  one  or  the  other  of  the  two  phases 
to  a  particular  clipper  amplifier  input  would  assure  that  the  corre- 
sponding shift  register  outputs  were  identical,  either  off  or  on. 
By  adjusting  the  number  of  clipper  amplifier  inputs  of  each  phase, 
the  matrix  output  for  any  beam  could  be  varied  from  0  to  +  16  units. 
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Detector  bias  was  adjusted  so  that  with  17  inputs  of  one  phase  and 
15  of  the  opposite  phase,  the  deteetor  output  amplitude  was  1/2  the 
difference  in  clock  pulse  amplitudes  as  viewed  by  oscilloscope  and  as 
shown  in  Fig,  10b.  Transient  spikes  as  indicated  in  the  figure  and 
detector  ineffieiency  prevented  the  achieving  of  exactly  aero  output 
voltage  for  the  oondition  of  inputs  divided  equally  as  to  phase. 

Time  variations  in  individual  detector  characteristics  were 
compensated  for  by  trimming  each  detector  bias  to  permit  alignment 
with  a  predetermined  test  pattern  as  shown  in  Fig.  10c.  This  was 
created  by  steering  the  test  signal  to  produce  maximum  response  on 
beams  1  and  J2. 

By  progressively  changing  the  inputs  from  half  of  each  phase 
(minimum  detector  output)  to  all  of  one  phase  (maximum  detector  out- 
put) complete  detector  characteristics  could  be  traced.  Many  trac- 
ings for  each  detector  were  required  ta  average  out  the  slight  dif- 
ferences in  individual  register  outputs  emphasized  by  the  artificial 
nature  of  the  single  frequency  input  but  averaged  automatically  in 
practice  by  the  random  nature  of  the  noise.  The  applicable  portion 
of  a  typical  curve  is  shown  in  Fig.  11.  After  one  set  of  character- 
istics had  been  determined  as  described,  a  duplicate  set  was  pro- 
duced. In  no  case  did  the  second  set  differ  by  more  than  1  millivolt 
from  the  first  set.  This  1  millivolt  error  is  considered  the  maxi- 
mum system  error  and  amounts  to  an  error  in  the  measurement  of  the 
output  of  any  beam  of  0.2  decibels. 
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5.5  The  Array 

Individual  array  elements  compricing  the  BaTiO^  receiving 
unit  and  the  FM  transmitter  canister  housing  were  calibrated  by  the 
U.  3.  Navy  Electronics  Laboratory  calibration  facility  at  S./eetwater 
Lake,  San  Diego,  California.  Measured  horizontal  and  vertical  re- 
ceiving patterns  were  completely  circular.  Since  the  wavelength  at 
the  nominal  array  center  frequency  of  1  kilocycle  per  second  is  5 
feet,  the  patterns  are  understandably  uniform  indicating  that  the 
individual  hydrophones  are  equivalent  to  a  point  hydrophone.  The 
cylindrical  shape  of  the  hydrophone  and  canister  does  not  alter  the 
isotropic  response. 

Although  the  response  of  the  assembled  array  of  individual 
elements  oould  not  be  measured,  acoustic  theory  and  '.Vorld  Var  II 
measurements ^  */  of  sonar  dome  structures  demonstrated  that  support 
members  of  cross-section  equal  to  or  less  than  1/15  ©£  a  wavelength 
have  no  effect  upon  the  beam  patterns  of  arrays  transmitting  through 
such  structures.  Since  the  ambient  noise  array  support  members  are 
1  inch  steel  angle  iron  approximately  1/6  inch  thick,  the  above  de- 
sign criterion  is  more  than  adequately  satisfied.  It  is  concluded, 
therefore,  that  the  array  may  be  considered  formed  of  ?2  exactly 
identical  omnidirectional  receiving  elements  unaffected  acoustically 
by  any  supporting  structure.  Accordingly,  it  is  considered  that  no 
errors  are  introduced  by  the  array. 
5.4  The  RF  System 

The  RF  system  was  examined  as  to  circuit  noise,  dynamic 
range  and  the  introduction  of  undesired  phase  shifts. 


31 


With  an  appropriate  capacitor  aubgtituted  for  the  hydrophone 
crystal,  electronic  circuit  noise  levels  averaged  20  decibels  below 
sea  state  zero  noise  levels. 

The  question  of  dynamic  range  primarily  concerned  channel 
cross-modulation.  Modulation  levels  were  adjusted  such  that  maximum 
noise  levels  as  measured  at  sea  produced  no  more  than  a  +  $0   kilo- 
cycle per  second  deviation  of  the  appropriate  channel  carrier.  This 
adjustment  permitted  a  100  kilocycle  per  second  separation  between 
channels  and  eliminated  all  possibility  of  cross-talk  between  chan- 
nels. 

Signal  phase  shifts  unintentionally  introduced  by  receiving 
channel  components  are  equivalent  to  eleotrical  time  delays.   Only 
differences  between  channel  phase  shifts  are  significant  since  equal 
phase  shifts  in  each  channel  are  equivalent  to  the  application  of  a 
fixed  time  delay  to  all  ohannelB  which  would  not  affect  beam  output, 
Channel  phase  shift  differences  will  not  seriously  affect  beam  out- 
put until  they  approximately  equal  the  sample  period  of  0,2  milli- 
seconds or  72  degrees.  A  large  relative  phase  shift  would  cause  the 
signal  to  appear  at  the  incorrect  shift  register  tap  and  would  dis- 
tort the  beam  psrttern.  To  measure  the  actual  shift  introduced  by 
the  system  from  hydrophone  to  shift  register  input,  conventional 
Lissajous  figure  osoillosoope  techniques  were  employed.  An  audio 
oscillator  output  was  coupled  by  suitable  attenuators  to  the  oscillo- 
scope horizontal  deflection  plates,  A  parallel  oscillator  output 
was  capacitatively  coupled  to  the  hydrophones  in  air  by  a  special 
cup  designed  to  fit  concentrically  over  the  array  hydrophone 
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©laments.  The  clipper  amplifier  output  of  the  particular  hydrophone 
channel  undergoing  measurement  was  fed  to  the  vertical  deflection 
plates  of  the  oscilloscope.  The  oscillator  frequency  was  adjusted 
for  a  straight  line  or  zero  phase  shift.  For  the  }2  channels  the 
frequency  of  zero  phase  shift  was  found  to  vary  between  +  JO   cycles 
per  second  centered  at  1.01  kilocycles  per  second.  A  second  measure- 
ment of  the  phase  shift  characteristic  of  an  individual  channel 
determined  by  observing  the  frequency  corresponding  to  the  easily 
read  Liesajous  points  of  180,  90  and  zero  degrees  permitted  conver- 
sion of  the  initial  zero  phase  shift  frequency  measurement  to  phase 
shift  at  zero  frequency.  From  thi6  latter  auxiliary  measurement  it 
was  established  that  0.35  degrees  shift  occurred  per  cycle  or  that 
the  greatest  measured  shift  was  4  10  degrees  from  1.01  kilocycles 
per  second,  a  value  well  within  the  sampling  period. 

It  is  concluded,  therefore,  that  no   significant  errors  were 
introduced  by  the  HF  system. 

^.5  Beam  Forming  System 

Errors  in  beam  output  introduced  by  the  beam  forming  system 
may  occur  due  to  four  possible  causes: 

a.  Inadequate  clipper  amplifier  dynamic  range. 

b.  Nonuniformity  of  clipper  amplifier  input  band-limiting 
filters. 

c.  Amplitude  variation  of  shift  register  outputs. 

d.  Nonuniformity  of  matrix  summing  resistors. 
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Clipper  amplifier  dynamic  range,  defined  as  that  range  where 
clipped  output  is  independent  of  input  amplitude,  was  measured  to  be 
70  decibels,  a  value  well  in  excess  of  the  range  of  ambient  noise 
levels  in  the  ocean.  No  error  is  attributed  to  this  source*  Ander- 
sons 5)  has  shown  that  errors  in  output  of  individual  elements  com- 
prising an  array  degrade  the  array  signal-to-noise  output  ratio  in 
accordance  with  the  relation: 

S/N  =  n  ( I  -  a2) 
where  n    is  the  number  of  array  elements  and  a   is  the  standard  devia- 
tion of  the  individual  element  errors.  The  standard  deviation  for 
items  (b),  (c)  and  (d)  above  when  measured  individually  and  then 
combined  was  less  than  10  per  cent.   Substitution  of  this  value  into 
the  above  expression  yields  a  0.05  decibels  loss,  a  negligible  con- 
tribution to  overall  system  error  from  beam  forming  components. 

5.6  Film  Recording 

As  noted  in  paragraph  2.3,  individual  beam  outputs  were 
photographically  recorded  continuously  throughout  each  run  by  l6mm 
camera.  Frame  spacing  was  geared  down  to  5*5  frames  per  second.  The 
developed  film  was  viewed  for  reading  on  a  Craig  16mm  Projecto  Edi- 
tor. The  displayed  frame  size  was  as  shown  in  Fig.  12  where  a  40 
millivolt  calibration  pedestal,  the  base  line  and  52  beam  outputs 
are  indicated.  An  accurate,  fine  line,  overlay  grid  of  A  millivolts 
per   division  was  drawn  to  large  scale  and  then  photographically  re- 
duced to  the  size  of  Fig.  12.  The  reading  procedure  first  involved 
adjusting  the  projected  frame  size  such  that  the  calibration  pulse 
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height  corresponded  exactly  to  AO  nillivolt9  on  the  viewing  grid. 
The  value  of  each  beam  output  cculd  then  be  read  quite  simply  to  the 
closest  r.illivolt.  Assuming  that  individual  component  errors  add 
randomly,  the  error  introduced  in  the  film  reading  process  is  con- 
sidered small  in  comparison  to  the  detector  calibration  error  of 
paragraph  5«2» 

5«7  System  Linearity 

A  consequence  of  digital  array  phasing  by  the  DIMUS  technique 
and  symmetry  in  the  array  beams  ^is  that  the  sum  of  the  output  powers 
of  the  J2  array  beams  is  a  constant.  Thus,  an  increase  in  level  in  a 
given  beam  will  result  in  a  corresponding  decrease  in  level  in  the 
remaining  beam  outputs,  A  very  strong  signal,  or  its  equivalent, 
good  correlation  in  a  given  direction,  can  result  in  beam  capture 
with  resultant  suppression  of  the  remaining  beam  outputs  and  attend- 
ant lo3s  of  information  as  to  the  relative  field  intensities  viewed 
by  several  beams.   It  is  necessary,  therefore,  to  determine  that 
maximum  output  level  on  a  given  beam  which  does  not  introduce  dis- 
tortion in  the  outputs  of  the  remaining  beams.  Toward  that  end  it 
is  assumed  that  a  signal-to^noise  ratio  at  the  input  of  an  individual 
hydrophone  of  1  to  10  represents  a  sufficiently  small  signal  such  as 
to  preclude  the  possibility  of  signal  capture.  The  output  power  gen- 
erated by  such  a  signal-to-noise  ratio  has  been  derived  in  Appendix 
II  and  the  results  compared  to  maximum  outputs  observed  during  mea- 
surements at  sea.   Since  no  measured  beam  outputs  exceeded  those 
computed,  it  is  concluded  that  the  system  behaved  in  a  linear  manner 
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throughout  the  experiment. 


IV.   TEST  PROCEDURE 
4.1   Narrative 

The  on-station  test  procedures  were  straightforward  barring 
the  occasional  equipment  component  failure.  Two  hours  prior  to 
commencement  of  test  the  equipment  was  energized  in  order  to  permit 
a  stable  temperature  to  be  reached.  The  detector  circuits,  for  ex- 
ample, proved  to  be  extremely  sensitive  to  temperature  changes 
greater  than  1  degree  and  required  thermocouple  control  of  the  two 
electronic  equipment  hut  cooling  fans.  Upon  reaching  the  test  area 
the  parachute  sea  anchor  was  rigged  and  the  ship  permitted  to  reach 
an  equilibrium  condition  headed  into  the  seas.  The  array  was  then 
swung  clear  of  the  ship  and  lowered  to  within  approximately  80  feet 
of  the  desired  operating  depth.  The  spar  buoy  and  float  system  was 
then  clamped  to  the  array  cable  and  the  system  lowered  until  the 
weight  of  the  array  was  supported  by  the  floats.  The  array  was  po- 
sitioned generally  between  100  and  1^0  feet  from  the  ship. 

Radar  bearings  and  ranges  to  any  contacts,  anemometer  read- 
ings and  fathometer  soundings  were  recorded  prior  to  turning  off  all 
ship  power  with  the  exception  of  the  ambient  noise  equipment  genera- 
tor. The  detectors  were  balanced  to  the  calibration  signal  as  des- 
cribed in  paragraph  5«2  and  an  approximate  4  minute  film  recording 
made  of  background  noise.  Array  compass  heading,  ship's  head,  swell 
direction,  bearing  to  the  array,  oscilloscope  and  calibration  pedes- 
tal and  base  line  settings  were  recorded. 

Normally  2  recordings  over  an  hour  interval  were  made  at 
each  depth.  At  the  conclusion  of  a  given  sequence  of  recordings  at 
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one  depth,  ship's  power  equipment  was  energized,  the  array  partially 
raised  to  permit  detaching  the  float  system  and  then  lowered  to  the 
new  depth.  The  procedure  as  described  above  was  then  repeated.  Ap- 
proximately 6   hours  were  required  to  complete  measurements  at  one 
station  of  5  array  depths. 

4.2  Test  Environment 

Tests  at  sea  were  planned  with  the  objective  of  maintaining 
a  constant  environment.  Since  ocean  surface  effects  as  noted  in 
Section  I  were  a  likely  source  of  non-uniform  ambient  noise  distri- 
butions, an  effort  was  made  to  control  other  possible  variables  such 
as  the  effects  of  land  masses,  surf,  shallow  water  and  varying  types 
of  ocean  bottom.  Although  it  would  have  been  desirable  with  the 
above  objective  in  view  to  conduet  all  tests  in  the  truly  deep  ocean 
beyond  the  continental  shelf,  only  one  such  operation  wag  undertaken 
owing  to  the  limitations  in  ship  speed  which  dictated  a  24  hour  one- 
way travel  time  to  such  areas.  Only  one  such  extended  cruise  was 
undertaken  and  accounted  for  all  sea  state  three,  2000  fathom  ocean 
depth  data  to  be  reported.  The  sea  state  one  plus  data  were  obtain- 
ed in  1000  fathom  water  in  an  area  which  shoaled  to  500  fathoms  in 
all  directions  within  approximately  10  miles.  The  sea  state  one- 
half  data  were  obtained  in  a  1000  fathom  area  similar  to  that 
studied  under  sea  state  one  plus. 

Sea  state  classifications  are  at  best  an  approximation  even 
when  derived  from  a  stable  platform.  Since  the  difficulty  of  main- 
taining an  objective  viewpoint  from  a  rolling  and  pitching  75  foot 
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test  platform  was  appreciated,  the  initially  recorded  sea  state 
values  of  two  and  four  were  reduced  to  one  plus  and  three  respec- 
tively following  examination  of  motion  piotures  taken  at  the  time  of 
recording.  All  observations  were  conducted  in  moderate  swell  con- 
ditions. The  environmental  conditions  are  summarised  as  follows: 

Sea  State    Ocean  Depth    Array  Depth    Type  Bottom 


5 

2000  fathoms 

150.  560, 
1000  feet 

Gray  mud 

H 

1000  fathoms 

--   150  feet 

Gray  mud 

1/2 

1000  fathoms 

1000  feet 

Gray  mud 

?.  MEASURED  DATA  ANALYSIS  AND  RESULTS 
5.1   Noise  Source  Origin 

Non-ieotropic  noise  distributions,  functions  of  both  azimuth 
and  elevation  angle,  were  measured  during  the  tests  described  in 
paragraph  4.2.  Observations  of  the  oharacter  of  the  noise  fluctua- 
tions in  the  upper  beams  of  the  array,  namely  the  vertical  upward 
beam  directed  normal  to  the  ocean  surface  and  the  next  family  of  5 
beams  at  58  degrees  from  the  vertioal,  indicate  that  one  of  the  sev- 
eral noise  generating  mechanisms  is  the  white  caps  of  breaking 
waves.  Near  the  ocean  surface  individual  beam  outputs  vary  consid- 
erably from  one  another  in  a  random  manner.  Momentary  increases  in 
level  in  a  given  beam  were  probably  due  to  the  arrival  of  energy 
from  an  individual  or  small  group  of  white  caps  or  wavelets.  As 
array  depth  increased,  beam-to-beam  differences  decreased.  The 
logical  explanation  lies  in  the  increased  area  of  ocean  surface  sub- 
tended by  the  beams  and  the  liklihood  that  the  beams  of  a  particular 
elevation  angle  view  more  nearly  the  same  distribution  of  white  caps. 
The  above  conclusion  was  substantiated  during  tests  designed  to  re- 
duce external  noise  interference  described  in  paragraph  2.6. 

Despite  efforts  to  sound  isolate  the  15  KW  power  generating 
equipment  during  early  phases  of  the  experiment,  interference  from 
the  test  ship  Invariably  was  observed  in  that  beam  most  nearly  di- 
rected at  the  ship.  When  installation  of  special  sound  isolated 
generating  equipment  failed  to  reduce  the  interference,  further  ex- 
amination revealed  that  the  test  ship  was  unwittingly  generating  the 
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equivalent  of  white  caps.  The  heavy  handling  boom,  rigged  to  port 
for  all  tests  created  a  pronounced  list.  The  list,  coupled  with  a  1 
to  2  foot  freeboard  amidships,  resulted  in  the  deck  edge  dipping  be- 
low the  waves  as  the  ship  rolled  even  in  the  lowest  sea  states.  A 
considerable  volume  of  water  entered  through  deck  edge  scuppers.  On 
the  return  roll  this  water  would  pour  from  the  scuppers  on  to  the 
ocean  surface  like  water  from  a  pitcher.   The  process  was  repeated 
with  each  roll.  The  interference  was  substantially  reduced  and  the 
source  identified  thereby  by  streamlining  the  entire  portside, 
rounding  over  projections  and  covering  scuppers  with  a  40  foot  long, 
4  foot  wide  strip  of  canvas. 

White  caps  are  probably  the  principal  contributor  to  back- 
ground noise.  However,  the  ocean  surface  by  some  unexplained  mech- 
anism still  contributes  to  the  noise  field  at  sea  states  above  zero 
even  when  white  caps  are  not  present.  This  fact  is  substantiated  by 
the  measured  distributions  reported  in  paragraph  6.% 

%2  Data  Analysis 

The  film  data  was  first  reduced  to  tabular  form  as  described 
in  paragraph  5»6.  Three  equally  spaced  4  minute  samples  recorded 
within  a  1  hour  interval  and  with  an  observation  every  14  seconds 
within  each  sample  resulted.  The  5  samples  for  each  array  position 
were  examined  for  fluctuations  within  each  sample  and  between  sam- 
ples. Fluctuations  within  samples  varied  from  less  than  +  0.4  deci- 
bels for  beams  directed  in  the  lower  hemisphere  to  ♦  1.0  decibels 
for  beams  directed  upward.  The  standard  deviation  of  sample 


42 


averages  was  smaller  generally  than  the  standard  deviation  of  obser- 
vations within  a  sample  by  the  factor  -jr    y2        where  n      is  the 
number  of  observations  comprising  a  sample.  Accordingly  three  sam- 
ples at  each  depth  were  averaged  and  treated  as  one. 

%3  A2lmuthal  Dependence 

The  data  were  first  analyzed  for  azimuthal  dependence  by  ex- 
amining the  outputa  of  families  of  beams  of  a  particular  elevation 
angle.  In  general,  some  azimuthal  dependence  was  apparent  princi- 
pally in  the  outputs  of  those  anray  beams  directed  in  the  upper 
hemisphere.  The  distributions  normally  were  elliptical  in  shape. 
A  representative  asimuthal  distribution  is  shown  in  Fig.  1%  It  is 
eonsidered  significant  that  the  maximum  beam  output  for  a  given  ele- 
vation angle  normally  aligned  itself  with  the  wave  fronts  of  the 
swells  while  mssrimum  outputs  usually  were  perpendicular  to  the 
swells.  This  fact  suggests  a  non-uniform  radiation  pattern  of  the 
white  cap  or  wavelet  source;  transmission  in  directions  normal  to 
the  wave  front  is  impeded  by  the  troughs,  with  lateral  transmission 
along  the  ridges  being  less  affected. 

5.4  Vertical  Dependence 

The  data  were  examined  for  n on- uniformity  in  the  vertical 
plane  by  averaging  the  outputs  of  all  beams  of  a  given  elevation 
angle,  thus  eliminating  azimuthal  dependence.  The  results  are  shown 
in  Fig.  14  which  presents  the  measured  field  as  a  function  of  eleva- 
tion angle  for  the  conditions  listed  in  paragraph  4.2.  Sea  state 
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FIGURE  13      AMBIENT  NOISE   AZIMUTHAL  DEPENDENCE 
FOR   VARIOUS    ELEVATION  ANGLES(0) 

SEA    STATE  3     ARRAY  DEPTH   130  FT. 
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FIGURE    14     VERTICAL    DISTRIBUTION    OF  AMBIENT  NOISE 
VS    SEA     STATE    AND    ARRAY    DEPTH 
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effects  are  apparent  with  maximum  to  minimum  differences  for  a  given 
condition  varying  from  less  than  1  decibel  for  sea  state  one-half  to 
6   decibels  for  the  highest  sea  state.  Less  significance  can  be  at- 
tributed to  differences  in  anisotropy  magnitude  due  the  variable  of 
array  depth  because  of  uncertainty  as  to  the  constancy  of  sea  state 
during  the  recording  period.  Approximately  6   hours  were  required  to 
record  data  at  the  three  depths  during  an  early  morning  period  of 
variable  winds.   It  is  not  unlikely  that  the  recorded  sea  state  of 
three  could  in  fact  have  fluctuated  between  sea  state  two  plus  and 
three  plus. 


VI.   TRtE  FIELD  RESTORATION 

6.1   True  Field  Indeterminacy 

In  the  preceding  paragraphs  the  analysis  has  been  confined 
to  the  measured  fields.   It  is  necessary  to  distinguish  between  the 
measured  distributions  and  the  true  noise  distributions.  The  output 
from  a  beam  whose  axis  is  oriented  in  a  fixed  direction  is  not  de- 
termined solely  by  the  true  field  intensity  in  the  oriented  direc- 
tion. The  beam  output  is  proportional  to  the  weighted  mean  of  the 
intensities  in  all  directions.  The  weighting  factor  is  the  beam 
pattern.  Mathematically  expressed,  the  output  of  a  beam  becomes: 

y^by~[  g\ob-  &,    f-AxWdn  (1) 

where     Y  m  beam  output  or  measured  field. 

G   ■  the  beam  response,  a  function  of  the  angular  difference 
between  axis  orientation,    9B  ,    <fig  ,   and  a  direction  of 
interest    <^\  <fi 
X  »  the  true  field  intensity  where,   for  simplicity,   aaimuthal 
symmetry  in  the  noise  distribution  is  assumed. 
Since  the  array  beam  patterns  are  not  pencil-like  with  aero  response 
in  all  but  the  desired  direction  but  are  broad  and  possess  a  side- 
lobe  structure  as  shown  in  Fig.  8,  the  effect  is  to  snear  the  fea- 
tures of  the  true  field. 

The  problem  of  reconstructing  the  true  from  the  measured 
distribution  is  common  to  the  field  of  astronomy  where  a  radio  an- 
tenna  is  used  to   survey  the  distribution  of  radio  brightness  over  the 
sky.     The  arguments  of  Braoewellv    '*'  who  first  considered  the  problem 
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in  detail  in  conjunction  with  aerial   smoothing  in  radio  astronomy  may- 
be adapted  to  the  underwater  acoustics  case  as  follows.     Upon  inte- 
grating oyer  the  variable,    0    ,   equation  (1)  becomes 

ri*B)=f      G{0B-  8)X\0)d&  (2) 

J-rr/z 

where  G    signifies  the  average  of  G   over  all  asimuth  for  a  given 
colatitude  and  X'  =    X  (^)sino?.  The  integral  on  the  right,  a  con- 
volution integral,  may  be  written 

Y{s)   =  G*X  (5) 

By  means  of  the  convolution  theorem,  equation  (5)  may  be  rewritten  in 
terms  of  the  Fourier  transforms  of  the  three  variables  as  follows 

?(s)    =  G(s)X(s)  (h) 

where  the  convolution  symbol  has  been  replaced  by  the  algebraic  pro- 
duet  and  the  three  terms  are  the  spectrums  of  the  corresponding  func- 
tions of  equation  (5)  and  5  '  is  the  number  of  cycles  per  unit  of 
angle,  &    .  From  (4)  it  is  seen  that  Y ( s )  =  0  whenever  G{s)-0 
and  thus  X  is)     is  indeterminate  under  such  conditions.  Obviously 
in  any  restoration  process  dependent  upon  measured  values  of  Y    only 
the  low  frequency  components  of  X   are  restored  to  their  full  value. 
The  resulting  values  of  X    have  been  called  the  "principal  solution" 
which  may  be  defined  as  that  solution  whose  spectrum  is  the  same  as 
the  spectrum  of  the  true  distribution  at  all  values  of  S  for  which 
G  is)   ^  0 .  The  results  of  a  restoration  proeess  to  be  described  in 
the  paragraphs  following,  are  shown  in  Fig.  1%  The  significant 
feature  of  Fig.  15  i«  that  the  transform  of  the  restored  field 
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FIGURE  15.  DISTRIBUTION    TRANSFORMS 
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approaches  zero  at  almost  the  same  rate  as  the  transform  of  the  mea- 
sured field.  This  is  an  important  condition  since  any  values  otX{s) 
for  s  greater  than  sk  t  the  value  of  s  where  G{s)=  0  ,  would  rep- 
resent fictitious  invisible  distributions  introduced  by  the  restora- 
tion process.  Since  X{sA)  ta  0  ,  the  addition  of  fictitious  high 
frequency  components  is  negligible  and  X  {0)  may  be  considered  an 
exact  principle  solution  throughout  this  study. 

6.2  The  Matrix  Approximation 

The  restoration  process  ^employed  can  be,  considered  in  two 
partsi   the  estimation  of  a  first  approximation,  and  an  iterative 
solution  deriving  from  such  first  approximation.  The  former,  the 
matrix  approximation,  is  novel  and  is  considered  in  thi3  paragraph. 
The  latter,  similar  to  methods  of  Bracewell,  is  described  in  para- 
graph 6.% 

The  matrix  approximation  is  generated  by  converting  the  con- 
volution integral  of  paragraph  6.1  to  a  finite  summation  of  terms  of 
number  equal  to  the  number  of  beams.  Each  term  in  the  summation  rep- 
resents a  contribution  from  a  particular  solid  angle  and  is  propor- 
tional to  the  average  of  the  true  field  subtended  by  the  solid  angle. 
Since  a  similar  expression  can  be  written  for  each  beam,  a  square 
matrix  results  relating  values  of  measured  field  to  the  unknown 
values  of  the  true  field.  By  inverting  the  matrix,  the  true  field 
values  can  be  obtained.  An  exact  matrix  can  be  written  only  for  the 
isotropic  case.  However,  the  matrix  so  derived  when  used  with  non- 
isotropic  distributions  yields  true  field  values  which  although 
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inexact,  are  closer  approximations  to  the  true  field  distributions 
than  can  be  arrived  at  by  intuition.  The  approximation  so  derived 
may  be  then  conveniently  operated  upon  by  iterative  techniques  to 
yield  the  principal  solution.  The  detail  a  of  the  derivation  of  the 
matrix  approximation  are  contained  in  Appendix  III. 

6.3  The  Iterative  Solution 

The  iterative  solution  is  a  method  of  restoring  the  true 
field  by  successive  substitutions.  The  first  approximation  pre- 
viously described  is  tested  by  numerically  completing  the  integra- 
tion of  equation  ( 1 )  of  paragraph  6.1,  and  comparing  the  results  with 
the  measured  values  of  beam  power,  Yl     .  From  experience  with  num- 
erous computations  of  this  nature,  it  has  been  determined  that  a 
suitable  second  approximation  is  obtained  by  correcting  the  first 
approximation  by  an  amount  equal  to  twice  the  difference  in  decibels 
between  the  measured  and  computed  values  of  Y/     •  Three  6uch  itera- 
tions were  usually  sufficient  to  reduce  the  difference  between  mea- 
sured and  computed  values  to  less  than  experimental  measurement 
error. 

Bracewell^ '    '   has  shown  that  a  sufficient  oondition  cf  con- 
vergence of  the  above  process  it  that  £(s)be  non-negative,  a  condi- 
tion which  is  satisfied  by  the  transforms  of  the  beam  patterns  of  the 
experimental  array. 

A  comparison  of  the  decibel  differences  between  a  given  pair 
of  beams  before  and  after  restoration  show  a  two  to  three-fold  in- 
crease in  the  decibel  difference.  It  is  reasoned  that  measurement 
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errors  will  be  similarly  increased  in  the  restoration  pr oeess.  Ac- 
cordingly, the  errors  in  the  true  field  principal  solution  are  es- 
timated to  be  no  greater  than  +0.6  decibels  except  in  those  regions 
where  abrupt  discontinuities  in  the  distributions  exist*  A  sample 
computation  of  the  complete  restoration  process  is  contained  in 
Appendix  17, 

The  resulting  true  field  distributions  for  the  five  environ- 
mental conditions  described  in  paragraph  4.2  are  shown  and  tabulated 
in  Pigs.  16  through  20.  The  existanee  of  a  marked  noise  null  in  the 
horizontal  direction  is  especially  useful  in  the  design  of  sonar 
arrays  whose  objective  is  to  detect  weak  signals  in  a  noise  back- 
ground. For  the  detection  of  horizontally  arriving  signals  in  a 
background  of  noise  arriving  from  all  directions,  good  vertical 
directivity  as  well  as  horizontal  directivity  is  advisable.  To 
illustrate!  consider  two  arrays  of  equal  horizontal  directivity  but 
differing  vertical  directivity.  Let  one  array  have  uniform  vertieal 
response  and  the  second  20  decibel  rejection  over  all  vertical  angles 
except  within  15  degrees  of  the  horizontal.  The  signal-to-noise 
ratio  of  the  latter  array  will  be  5*5  and  8  decibels  better  than  the 
former  when  both  are  positioned  in  an  isotropic  noise  field  and  the 
field  of  Fig.  16  respectively. 

As  previously  noted,  variations  in  the  noise  field  distri- 
bution for  differing  array  depths  could  not  be  established  because 
of  uncertainty  as  to  the  constancy  of  sea  state  during  the  recording 
period.  The  increase  in  anisotropy  with  increasing  sea  state,  how- 
ever, is  obvious. 
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FIGURE  16    TRUE  FIELD  VERTICAL  DISTRIBUTION 
SEA    STATE  3     ARRAY  DEPTH   1000' 
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FIGURE  17    TRUE  FIELD  VERTICAL  DISTRIBUTION 
SEA    STATE  3     ARRAY  DEPTH  560  FT 
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FIGURE  18     TRUE  FIELD    VERTICAL   DISTRIBUTION 
SEA    STATE  3    ARRAY    DEPTH    150    FT. 
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FIGURE  19     TRUE   FIELD   VERTICAL  DISTRIBUTION 


SEA   STATE   I     ARRAY   DEPTH    130   FT. 
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VII.  AN  OCEAN  MODEL 

7*1  Ocean  Model  Development 

From  the  nature  of  the  distributions  shewn  in  Figs,  16 
through  20  and  consistent  therewith,  the  following  simple  picture  of 
the  mechanism  of  noise  generation  and  redistribution  in  the  ocean 
has  been  constructed.  Noise  in  the  deep  ocean  is  presumed  to  be  a 
superposition  of  two  fields.  The  first  field  is  that  which  exists  at 
sea  state  zero  prior  to  the  development  of  wave  motion  and/or  white 
caps.  Its  distribution  is  isotropic  and  its  origin  was  not  estab- 
lished during  this  experiment.  A  possible  explanation  of  the  origin 
of  low  sea  state  ambient  noise  is  considered  in  paragraph  8.1.  The 
second  source  of  noise  is  generated  by  wave  motion  and  white  caps  at 
the  ocean  surface.  The  equivalent  surface  radiator  is  highly  direc- 
tional in  the  vertical  plane  when  asimuthal  effects  are  averaged. 
The  directionality  can  possibly  be  attributed  to  the  ocean  waves  and 
swells  functioning  in  a  manner  analogous  to  angular  focusing  reflec- 
tors. In  addition,  in  azimuth  the  ocean  swells  impede  transmission 
at  shallow  transmission  angles  in  directions  normal  to  the  wave 
fronts  of  the  swells. 

It  is  further  postulated  that  the  noise  intensity  at  a  given 
point  in  the  ocean  due  to  the  equivalent  surfaoe  radiator  arrives  by 
both  a  direct  ray  path  and  single  bounce  path  experiencing  specular 
reflection  at  the  ocean  bottom.  The  intensity  at  the  surface  souroe 
experiences  inverse  square  spreading  loss  and  absorptive  attenuation 
of  the  form  10  10  where  a  is  the  absorption  coefficient 
( decibel s/kiloyard)  and  R   is  ray  path  length  in  kiloyards. 
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Referring  to  the  model  of  Fig,  21   and  disregarding  for  the  moment 
the  lov  sea  state  isotropic  contribution  te  the  total  noise  field, 
the  abore  postulation  of  an  ocean  model  may  be  restated  in  mathe- 
matical form  as  follows: 

-aRp 
X(9)df2-     S{*\dA     |0     '0  (1) 

vhere    X(9)dS2     «=  noise  intensity  per  incremental   solid  angle 

at  the  receiver  in  the  upper  hemisphere 
S{9)dA     =  noise  intensity  at  one  yard  per  incremental 
radiating  surface. 

Since 

d  A  cos#  (2) 


L 

equation  (1 ) 

may  be  rewri 

tten 

as 

X(9)dS2   =    ^^ 
cos/9 

dSl 

10 

-aRp 
10 

or,  finally, 

aH 

sec^ 

S(<?)  = 

X{9)  cos  S> 

10 

10 

(5) 


0  <  9  <    7T/z  (4) 

A  similar  expression  for  the  intensity  at  the  receiver  in  the  lower 
hemisphere  due  to  bottom  reflection  may  be  written 

-a{Rm+RR2) 
XiSL)dfi*   S{9)dA«{9)   |Q  ,o  () 

where  a    (  9)   i«  the  intensity  bottom  reflection  coefficient  and 
i  s  the  suppl  ement  of  9R  . 
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FIGURE    21         OCEAN  MODEL 
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Since 


dn-_      dA  cos  e  {6) 


equation  (5)  may  be  written  as 

a(ZD-H) 


sec  0 


a{6>)-_    X±lACos0\O  l0  O<0<    rr/Z  (7) 

S(0) 

tt/2  <  <9r  <   rr 


7.2  Deduced  5 (0)   and  a  (0) 

Employing  the  distributions  of  Figs*  16  through  19  and  equa- 
tions (4)  and  (7)  above,  values  of  S  {0)  &n&  a   (#)were  computed.  It 
was  necessary  to  estimate  the  contribution  of  the  isotropie  xero  sea 
state  noise  field  and  to  subtract  that  amount  from  the  total  field, 
X  (0)  ,  prior  to  proceeding  with  the  computation.  Several  alterna- 
tive methods  of  estimating  the  zero  sea  state  contribution,  none  of 
which  can  be  considered  exact,  are  possible.  For  example,  the  level 
of  the  null  occurring  at  an  elevation  angle  of  approximately  100  de- 
grees in  all  distributions  may  be  considered  due  to  the  low  sea 
state  isotropic  field.  Errors  in  this  assumption  are  likely  due  to 
the  difficulty  in  establishing  the  true  minimum.  The  most  probable 
point  of  error  in  the  restored  distribution  is  at  the  minimum  due  to 
the  difficulty  inherent  in  resolving  a  sharp  discontinuity  with  a 
broad  beam  pattern  array.  In  addition,  contrary  to  the  very  simple 
straight  ray  assvanptions  of  the  ocean  model,  curved  path  arrivals 
direct  from  the  surface  or  by  convergent  zone  paths^  '  dus  to  verti- 
cal velocity  gradients  probably  contribute  to  the  null  level.  For 
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example,  Bell^  ''  ha8  considered  the  deep  ocean  ease  with  negative 
gradient  to  deep  depth.  He  reports  bottom  reflection  losses  for 
grazing  angles  of  incidence  which  are  approximately  6   decibels  less 
due  to  the  refractive  effect  of  the  negative  gradient  than  would  be 
predicted  on  the  basis  of  straight  ray  path  assumptions. 

The  best  estimate  of  the  isotropic  contribution  utilized  the 
Knudsen^  '  sea  state/noise  level  curves.  If  the  restored  distribu- 
tion is  averaged  over  4  7r  S2   ,  the  relative  level  derived  from  an 
omnidirectional  hydrophone  can  be  obtained.  Subtracting  from  this 
level  the  decibel  difference  between  the  Knudsen  noise  level  for  the 
sea  state  under  consideration  and  sea  state  zero,  the  level  of  the 
isotropic  contribution  can  be  established.  The  weakness  in  the  above 
approach  rests  in  its  dependence  on  subjective  estimates  of  sea  state 
level.  The  error  introduced  by  this  technique  decreases  with  in- 
creasing sea  state  since  the  percentage  contribution  to  the  total 
field  of  zero  sea  state  isotropic  noise  becomes  smaller  under  such 
conditions. 

Figure  22  is  a  plot  of  Si  8),   the  surface  source  radiating 
pattern,  deduced  by  the  above  procedures  from  the  restored  fields  of 
Pigs,  16  through  20,  A  reasonable  approximation  to5(<^)is  the 

functioni 

90-<? 
S{0)  =  {cos  8)      l0  0   <     8  <   70° 

S{8)  =    cos2  8  70°  <   $  <  90° 

which  is  shown  superimposed  on  the  experimental  data. 
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FIGURE  22      SURFACE    SOURCE    RADIATION    PATTERN 
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The  bottom  reflection  coefficient  (reflected  pressure/ 
incident  pressure)  derived  for  conditions  corresponding  to  Fig.  22 
are  plotted  in  Fig.  23  together  with  experimental  data  and  theoreti- 
cal plots  from  MacKenzie^   '  and  Marsh'  ' '.  The  fit  is  considered 
good  in  view  of  the  indirect  manner  in  which  the  data  were  obtained 
and  the  several  assumptions  involved  in  the  analysis.  The  results 
are  considered  to  substantiate  the  simple  ocean  model  proposed. 
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VIII.  SUMMARY 

8.1   Conclusions  and  Discussion 

Prior  to  this  experiment,  as  noted  in  Chapter  I,  no  experi- 
mental data  existed  to  describe  the  directional  distrlbittion  of  deep 
ocean  ambient  noise.  Accordingly,  stated  objectives  of  this  study 
were  to  investigate  anisotropy  of  ocean  ambient  noise,  to  isolate 
vertical  and  azimuthal  dependencies,  to  establish  the  identity  of 
sources  contributing  to  non-uniform  distributions  and  to  infer  the 
radiation  patterns  of  such  sources.  All  the  above  research  objec- 
tives have  been  realized  by  this  experiment.  Deep  ocean  ambient 
noise  is  anisotropic  in  the  frequency  band  750  to  1500  cycles  per 
second  and  at  sea  states  above  2ero.  Vertical  and  asimuthal  depen- 
dencies have  been  isolated.  The  surface  of  the  ocean  ha6  been 
identified  as  a  major  contributor  to  non-uniform  distributions  and 
the  radiation  patterns  of  suoh  surface  sources  have  been  deduced. 

Vertical  anisotropy  arises  from  the  superposition  of  a 
noise  field  due  to  the  radiating  ocean  surface  upon  an  isotropic 
field  present  in  the  absence  of  surface  radiation.  The  shape  of  the 
vertical  distribution,  Figs.  16  through  20,  may  be  compared  to  the 
form  of  an  arrowhead.  The  arrow  tip  is  directed  normal  to  the  ocean 
surface  and  represents  maximum  intensity  arriving  from  the  surface 
direction.  The  arrowhead  tail  indioates  noise  originating  at  the 
ocean  surface  which  has  been  specularly  reflected  from  the  ocean 
floor.  A  dip  in  the  arrowhead  tail  in  a  direction  normal  to  the 
ocean  bottom  is  evidence  of  a  lower  value  of  bottom  reflection  co- 
efficient for  normally  incident  sound  than  for  other  angles.  A 
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narked  null  at  the  waist  of  the  arrowhead,  approximately  contained 
within  +  15  degrees  of  the  horizontal,  exists  in  the  vertical  dis- 
tribution. The  magnitude  of  the  null  varies  between  2  and  12  deci- 
bels below  the  maximum  intensity  normal  to  the  surface  for  sea  states 
one-half  to  three  respectively.  The  null  is  due  to  two  effects. 
First,  contributions  from  the  surface  source  to  the  horizontal  or 
near  horizontal  distribution  arrive  from  greater  distances  than 
contributions  to  the  off -horizontal  distributions  and  have  exper- 
ienced, therefore,  greater  attenuation  from  spreading  and  absorp- 
tion. Second,  the  level  of  radiation  from  the  surface  source  in 
horizontal  directions  is  more  than  20  decibels  below  the  level  of 
radiation  in  the  downward  directions. 

The  non-uniform  noise  distribution  with  a  minimum  in  the 
horizontal  direction  ia  of  practical  importance  to  the  design  of 
3onar  receiving  arrays.  The  importance  of  vertical  directivity  to 
the  detection  of  horizontally  arriving  signals  ia  apparent. 

The  origin  of  the  isotropic  zero  sea  state  noise  field  was 
not  established  by  this  experiment.  Of  the  various  sources  con- 
sidered by  Preach'  '  distant  surface  motion  is  considered  most 
likely.  Although  the  aea  surface  may  be  flat  at  the  point  of  mea- 
surement, the  noise  produced  by  distant  surface  disturbances  can 
propagate  very  long  distances  with  negligible  attenuation  due  to 
absorption  in  the  750  to  1J500  cycle  per  second  frequency  band  ob- 
served. The  scattering  and  reflection  processes  experienced  by  the 

sound  in  reaching  the  point  of  observation  could  destroy  its  origi- 
nal character  and  an  isotropio  field  result. 
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Ocean  swells  are  probably  responsible  for  azimuths!  aniso- 
tropy  by  impeding  transmission  of  surface  generated  noise  in  direc- 
tions nomal  to  the  wave  front  of  the  swells. 

The  intensity  at  a  receiver  due  to  a  unit  incremental  sur- 
face radiator  of  Intensity, S (*?),  at  one  yard,  can  be  considered  as 
the  sun  of  intensities  due  to  direct  transmission  find  specular  hot- 
ton  reflection  with  attenuation  due  to  spherical  spreading  and  ab- 
sorption occurring  over  the  separation  between  radiator  and  receiver. 
The  quantity,  S{&),   is  markedly  directional  in  the  downward  di- 
rection and  can  be  represented  approximately  by  the  expression: 

96--C 

(cos  0)  0  <    0  <    70° 

cos2^  70°  <  9<  90° 

The  above  ocean  model  is  not  unique.  It  does,  hov/ever,  con- 
form to  the  measured  data  as  evidenced  by  the  reasonable  agreement 
between  bottom  reflection  coefficients  determined  from  this  experi- 
ment and  values  reported  by  other  investigators.  Further  additions 
to  the  model,  as  for  example  the  consideration  of  bottom/surface 
multiple  reflections,  are  obvious.  Such  refinements  would  presume, 
however,  a  knowledge  of  bottom  topography  and  array  position  rela- 
tive to  the  bottom  greater  than  existed  during  the  experiment. 

Although  the  data  upon  which  the  ocean  model  is  based  were 
obtained  in  ocean  depths  greater  than  1000  fathoms,  there  is  no 
known  feature  of  the  model  which  would  preclude  its  application  to 
shallower  depths  of  perhaps  100  fathoms.  Its  application  to  even 
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shallower  depths  will  depend  upon  the  presently  unknown  effect  of 
water  depth  upon  the  nature  of  the  surface  source  radiation  pattern. 
On  the  basis  of  measurements  conducted  over  a  one  hour  per- 
iod under  constant  conditions  of  array  orientation  and  sea  state, 
the  ambient  noise  field  can  be  considered  stationary.  A  single, 
trainable  and  tiltable  directional  receiver  could  be  employed, 
therefore,  to  duplicate  the  performance  of  the  52  element,  multi- 
beam  array.  However,  to  achieve  beam  patterns  equivalent  to  those 
of  the  52  element  array  with  a  single  directional  unit,  focusing 
reflectors  of  dimensions  comparable  to  the  32  element  array  would 
be  required.  A  very  difficult  problem  of  mechanical  training  and 
elevating  the  reflecting  unit  would  be  substituted  for  the  compara- 
tively simple  telemetering  problem  associated  with  the  J2  element 
array. 

8.2  Suggestions  for  Further  Experiment 

The  results  contained  herein  represent  only  a  first  step  in 
the  area  of  ambient  noise  directional  distribution  investigation. 
Objectives  for  future  experiment  should  include t 

a.  Shallow  water  investigation. 

b.  Extension  to  higher  sea  states. 

e.  Extension  to  array  depths  of  at  least  5000  feet. 

d.  Extension  of  bandwidth  covered  from  less  than  A00 
cycles  per  second  to  5  kilocycles  per  second. 

e.  Investigation  of  the  effect  of  land  masses  such  as 
islands  and  the  continental  shelf  and  the  influence  of  surf. 
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An  additional  area  of  investigation  1b  to  substantiate  the- 
oretically and  experimentally  the  radiation  pettern  deduced  from  this 
experiment  of  a  white  cap  situated  on  the  crest  of  a  swell. 
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APPENDIX  I 

QUANTITATIVE  SYSTEM  ANALYSIS 

The  objective  of  this  section  is  to  derive  an  expreaaion  for 
the  voltage  appearing  at  the  output  of  a  given  beam  detector  for  a 
random  noiae  input. 

Oonaider  first  an  N  element  array  of  omnidirectional  elements 
situated  in  an  isotropic  noise  field  of  purely  continuous  spectrum. 
If  the  hydrophone  electrical  outputs  generated  by  the  noise  field  are 
summed,  with  or  without  the  introduction  of  an  electrical  delay,  and 
if  the  hydrophones  are  sufficiently  widely  spaced  such  that  cross- 
correlations  between  elements  are  zero,  the  array  output  following 

1/2 

square  law  detection  will  be  proportional  to  N        .If  the  cross- 


correlations  are  not  equal  to  zero,   the  output  will  be  proportional 

to    N  l/2   (  I  +  e  )  '/2      where 

.       /v(/v-D 

e  =   -rr  ZZ  p    (0) 
N     t*j      ,J 

and  p         is  the  normalised  cross-correlation  function  derived  from 
ij 

the  input  noise  spectrum.  The  value  of  e    may  also  be  determined 
from  the  directivity  ratio.  Utilizing  the  beam  patterns  of  paragraph 
2.6,  the  directivity  ratio  is  derived  by  numerical  integration  from 
the  expression 

4  w 


f      G[&,</>)dti 
where  G  {&,  ^)is  the  beam  pattern  of  paragraph  2.6,   and  6"  (0,0)  =     I 
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The  quantity,  "  £ "  ,  may  be  determined  from  the  expreesion 
/V/l  f  e   r  I  / D.  R.  .  For  inner  and  outer  beams  "  e"   has  been  com- 
puted to  be  0.25  and  0.02  respectively.  Lawson^1"^  has  shown  that 
the  correlation  functions  of  undipped  and  clipped  signals  are  re- 
lated by  p  =2/V  arcsin  P  •  For  omall  values  of  correlation 

'clipped  '^undipped 

as  is  the  case  under  consideration,  the  arc  sine  of  the  argument  is 
equal  to  the  argument  and  the  correlation  function  of  the  clipped 
signal  reduces  to  2  /V  times  the  undipped  correlation  function. 
Sinee  the  actual  detector  employed  was  a  linear  rectifier  rather  than 


a  square  law  device,  an  added  factor  K    -   0.8  is  required.  K--J2/t 
is  the  ratio  of  the  mean  absolute  value  to  the  root-mean- square 
value  in  a  normal  distribution.  A  further  factor  of  0.5  is  intro- 
duced because  the  scale  of  the  abscissa  of  Fig.  11  is  ^6   for  maximum 
detector  output  rather  than  J2.  The  detector  output  voltage  is 
therefore  proportional  to 

Substitution  of  the  values  N   »  52  and  e    -0.25  into  the  above  ex- 
pressions produces  the  value  2.44.  This  quantity  when  entered  as 
the  abscissa  of  Fig.  11  indicates  an  output  value  of  52  millivolts 
for  detector  No.  1,  an  inner  beam. 

The  computed  value  of  52  millivolts  was  experimentally  con- 
firmed from  data  recorded  at  sea.  When  the  outputs  of  all  J2  beams 
are  averaged,  the  resulting  value  is  equivalent  to  an  individual 
beam  output  to  be  expected  from  an  isotropic,  random  noise  field. 
The  averages  of  several  sets  of  data  equalled  approximately  52  milli- 
volts. 


APPENDIX  II 
SYSTEM  LINEARITY 

( 11  ) 
Rudniekv   '  includes  the  neoeasary  expressions  to  compute 

beam  output  power  a  a  a  function  of  array  geometry  and  input  signal- 
to-noise  ratio  at  an  individual  hydrophone.  Output  power  may  be 
expressed  as 

c    (0) 

(1) 


p   r  N    N  N     N 


l-\     J-\  l-\    /-   I 


with  the  summations  to  be  completed  over  the  N   array  elements.  The 
quantity,  M lf         (0)  ,  is  the  combined  autocorrelation  function  of 
signal  and  noise.  For  purposes  of  this  analysis,  the  spectra  of 
signal-to-noise  may  be  considered  identical  and  are  assigned  dif- 
ferent nomenclature  only  to  differentiate  between  that  noise  which  is 
within  the  main  lobe  of  the  directed  beam  and  that  which  is  outside 
the  lobe.   Implicit  in  the  last  statement  is  the  assumption  that 
noise  is  completely  uncorrelated  with  respect  to  angle  of  arrival. 
C      (J    (0)  is  the  autocorrelation  function  of  the  speotrum  introduced 
by  clipping  and  is  defined  by  Rudnick'   '  as 


2_ 

77- 


Sirf  u     (r) 


_2_ 

7T 


(2) 


ju      (  0 )     may  be  expanded  to 


ij 


M   .(0)   '-    -?—    fi     (0)  +    — ?—  a    (0) 
^ 'J  n-rs      '  n  +  s      'J 
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(5) 


75 


where  p    and  a   are  noiee  and  signal  autocorrelation  functions  re- 
spectirely.   er^    (0)  =  I   since  by  definition  it  is  the  signal  to 

N     V 

which  the  beam  is  steered.  Z  Z      c,,  (0)may  also  be  separated  into 


/--i  j 


N      N  N  N{N-\) 

Z  Z   c    (0)  =  Z  c,  (0)+  Z Z      c.i  (0)  (k) 

t=\    )-\       "■  /--I  i* j  '  v  ' 

It  follows  from  (2)  that      C;i    -  rr  /2  -|  •     In  addition,   Rudniek  has 
shown  that 

N    N 

ZZ   c.  (0) 

is  small  and  may  be  neglected. 

Substituting  the  several   expressions  described  above  into 
(1)  and  reducing  terms,   the  output  voltage  of  the  linear  half-wave 
rectifier  can  be  shown  to  be  proportional  to 


I/2 


^  2  rr'  L         sr     /?+5  rr     /7+ S    wv      '  'J  ^' 


l/2 


If  a  small  signal-to-noise  ratio  input  of  5 / n     equal  to  0.1 
is  assumed  and  with  N   equal  to  52  and  "<?"  equal  to  0.25,  equation 
(5)  yields  the  value  5*87.  Entering  the  detector  characteristic  of 
Fig.  11  with  the  above  value,  an  output  voltage  of  90  millivolts  is 
obtained.  No  recorded  data  was  used  in  the  analysis  where  the  out- 
put of  any  beam  exceeded  90  millivolts.  It  is  concluded,  therefore, 
that  the  system  was  operated  throughout  its  linear  region  during  the 
tests. 


APPENDIX  III 

DERIVATION  OF  THE  MATRIX  APPROXIMATION 

Considering  the  case  of  azimuthal  symmetry,  the  output  of 
any  vertical  beam  may  be  written » 

K  =  [  x{#)  g,  {&,<£)  dn  (1) 

where  Yf      «  output  of  / ,h     beam 

X     »  continuous  true  field  distribution  (power/unit 

receiving  area/untt  solid  angle  subtended  by  the 
source) 
G     a   power  response  of  / ,h   beam 
If  the  true  field  is  assumed  to  be  isotropic,  then 

r.  --  xn  f  gi  {&,<p)da  (2) 

The  direct! rity  ratio  of  an  array  is  defined  asi 

DR.    =  -&*- -. (5> 

c,  {#,  #)  f  dn 

Combining  equations  (2)  and  (5) 

Y/    =    Xn  G  tyrf)  OR  Arr  (4) 

where     G  ( 9  ,  <f>  )   equals  unity  and  0.853  for  outer  and  inner  beams 
respectively  as  defined  by  paragraph  2.6.     Setting    4  w  X^    -    I 
which  has  the  physical  interpretation  of  total  power  from  all  direc- 
tions normalised  to  unity,   then 


7* 
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D.R.    (outer  beam) 

Y,   =  (5) 

0.855  D.R.    (inner  beam) 

Rewriting  equation  (1)  at  a  sum  of  eight  terms,   with  each  term  rep- 
resenting the  contribution  from  a  spherical  zone  formed  by  two 
planes  parallel  to  the  equator  and  each  zone  containing  one  of  the 
eight  vertical  beams,   then: 


y,    ''  Z  *    J      s\r\0c/0  f    Gl{8y<fi)d<£ 


Jz\ 


'zone 


(6) 


The  true  field,  X  j   ,  may  be  non-isotropic,   although  X    is  a  constant 
within  a  particular  zone.     The  quantity    G  /  (&,0)  may  next  be  aver- 
aged  over  all  azimuth. 

J--8  M&/2  

y,    =  Z    X,     I        Zrr    G[#)     SM0d9  (7) 

/  =  .         '    % 


where    &  -^  and    &J2  are  the  limits  of  the  spherical   zones.     There 
exists  between  the  limits  of  eaeh  zone  a  particular  value  of @z@^\ 
which  permits  equation  (7)  to  be  transformed  to  the  following 

y,   *Z     Xj2wGt{0M  )    l^9d8  (8) 

/=  i  &)  i 

Equation  (8)  may  be  further  simplified  by  replacing  the  integral  by 

a  constant,  k     ,  for  a  particular  zone  yielding 

J-Q 

Y,    -  ZrG,   (*■„»)  Z      k      Xj  (9) 

l-\ 
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or  more  briefly 

y,   -"  atJ     Xj  i,  j  =   |-8  0°) 

Equation  (10)  relates  the  measured  field  values,  Y/  ,  to  a  set  of 
eight  simultaneous  equations  in  eight  unknowns  which  are  the  true 
field  ralues,  Xj   • 

The  a ,j     matrix  which  satisfies  equations  (5)  &n&   (9)  was 
determined  by  trial  and  error.  The  matrix  was  then  inverted  by 
digital  computer  to  produce  a  set  of  equations  of  the  form 

X  ~-   a'n      Y, 

The  coefficients  of  the  inverted  matrix  are  listed  below.  By  mul- 
tiplying each  element  of  row  one  by  Yi    ,   row  two  by  K,  ,  etc.,  the 
sum  of  all  oolumn  one  elements  yields  X x   ,  the  column  two  elements 
X z,   etc.  The  values  of  X   so  derived  are  the  first  approximation 
to  the  true  field  in  the  restoration  process. 

Table  of  Coefficients  a 

i)  1  2  5  4  5  6  7  8 

1  67.855  -52.557     27.269  -19.912     10.678  -19.522       5.122  -8.952 

2  -1.677  89.527  -67.679  25.059  -19.545  16.75*  -5.156  -1.845 
5       -O.O56     -9.622     92.556  -59.901      20.905  -52.688       0.808  -0.282 

4  -1.102  2.505  -52.047  106.681   -58.812     20.422     -5.628  -0.799 

5  -0.799  -5.628     20.422  -58. 812  106.681   -52.047       2.505  -1.102 

6  -0.282  0.808  -52.688     20.905  -59.901     92.556     -9.622  -O.O56 

7  -1.845  -5.158     16.754-19.545     25.059-67.679     89.527-1.677 

8  -8.952  5.122  -19.522     10.678  -19.912     27.269  -52.557  67.855 


APPENDIX  IV 

TRUE  FIELD  RESTORATION  SAMPLE  COMPUTATION 

In  order  to  illustrate  the  restoration  process,  the  proce- 
dures followed  in  deriving  Fig.  18  will  be  described. 

Table  I  ie  a  random  selection  of  recorded  data  typical  of 
that  which  results  following  the  film  reading  procedures  of  para- 
graph J.6.  The  values  shown  are  the  individual  beam  detector  out- 
puts in  millivolts  following  one  second  integration  and  at  14  second 
intervals.  Beam  numbering  is  ruch  that  beam  J2  is  directed  at  the 
surface,  beams  27-J1  at  58  degrees  from  the  vertical,  beams  22-26 
at  6b   degrees  from  the  vertical,  beams  17-21  at  78.5  degrees  from 
the  vertical,  etc.,  to  beam  1  directed  at  the  bottom  of  the  ocean. 
The  RMS  values  for  each  beam  for  each  run  are  computed  and  tabulated 
in  colurwi  2  (Aver)  of  Table  II,  The  detector  correction  is  applied 
(Corr  Aver)  and  the  RMS  output  for  beams  of  a  given  colatitude  com- 
puted (Vert  Aver).  The  average  of  all  beams  ie  similarly  determined. 

Power  ratios  for  each  colatitude  with  respect  to  beam  52  and 
with  respect  to  the  average  of  all  beams  are  computed.  The  Yf 
column  is  determined  by  multiplying  the  appropriate  expression  of 
equation  (5)  of  Appendix  III  by  the  power  ratios  of  column  8  of 
Table  II.  The  values  of  Y    are  then  used  in  conjunction  with  the 
matrix  of  Appendix  III  to  obtain  the  first  approximation  to  the  X  t   . 

In  Fig.  24  the  measured  values  of  Y/    ,  expressed  as  deoibel 
ratios  with  respect  to  beam  52,  are  plotted  as  dots.  A  faired 
curve  connecting  the  first  approximation  values  of  X ,   is  indicated 
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TABLE  I 
Beam 

1  51  50  51  50  51  51  51  52  50  51  51  50  50  50  51  50  51  51  50 

2  49  49  49  48  48  48  49  *9  *8  *9  48  48  49  49  48  46  47  49  49 
5  47  46  46  46  47  46  47  48  47  45  46  47  47  46  46  46  46  47  47 

4  49  49  48  48  49  48  48  49  49  48  49  49  49  49  49  48  49  49  49 

5  50  50  49  49  50  49  50  49  49  49  49  49  50  50  49  49  50  51  49 

6  47  47  47  47  47  46  47  46  45  47  47  46  46  47  46  46  47  46  47 

7  49  50  49  49  49  48  49  50  49  54  51  50  50  49  50  49  50  50  50 

8  50  51  50  50  49  50  50  41  50  50  50  50  51  49  50  50  51  50  50 

9  48  49  49  48  49  49  50  49  48  49  49  49  49  48  48  47  40  50  48 

10  49  52  50  50  50  50  51  52  50  50  50  51  50  50  50  51  51  52  50 

11  46  48  47  46  47  46  48  47  48  46  47  47  48  48  47  47  48  48  46 

12  42  44  45  44  45  44  45  45  45  45  44  44  45  44  44  45  44  45  45 
15  42  44  45  45  45  42  44  45  44  45  44  45  45  4}  45  45  44  45  44 

14  41  45  45  45  45  42  44  45  45  41  42  4}  42  42  42  42  4j  44  45 

15  46  49  49  49  49  48  49  50  48  47  49  49  49  49  48  47  48  49  47 

16  50  55  54  54  54  51  55  54  53  50  5^  51  52  52  52  52  51  52  52 

17  47  50  49  49  49  49  49  49  48  56  52  50  50  49  50  49  49  50  49 

18  54  58  51  56  56  54  57  56  56  56  57  56  56  55  55  56  56  56  54 

19  48  51  51  50  50  50  51  50  50  49  49  50  51  50  A9  48  50  52  49 

20  50  52  51  50  51  50  51  5^  52  51  50  50  49  50  51  51  50  50  51 

21  54  56  51  55  57  56  55  55  55  52  5*  5*  55  55  55  52  55  5*  51 

22  50  54  51  51  51  50  52  55  5^  52  52  55  52  52  52  50  52  52  51 
25  64  64  64  65  64  65  65  65  64  61  66  65  66  64  60  65  65  64  65 
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TABLE  I  (COOT.) 
Beam 

24  62  62  62  6^   61  61  62  6^   60   59  61  62  62  60  62   59  60  63  59 

25  51  53  50  51  50  50  50  52  51  50  52  51  52  51  51  51  51  51  51 

26  57  59  58  58  56  54  56  59  57  56   58  57  57  58  56  57  57  60  57 

27  60  59  59  60  58  58  59  58  58  56  58  58  59  59  58  57  58  60  57 

28  60  59  60  61  60  59  59  59  61  58  60  61  60  58  58  57  53  62  58 

29  66  62   62  65  62  6k  65  64  63  62  65  66  62  62   64  60  61  65  63 
JO  69  74  67  65  64  64  66   65  61  64  66  69  6l  68   63  64  70  63  63 

31  64  63  62  64  61  64  62  66  62  60  62   64  61  65  62  62  61  61  59 

32  77  78  76  74  74  76  75  77  73  74  75  76  72  74  75  76  77  75  72 


so 


TABLE  II 


Corr  Vert  Elev  Ratio  Ratio   Ratio     Ratio 
Bean  Aver  Aver  Aver*  Angle  re  52  re  av.  Re  32(db)  Re  av.(db)  Yt 


1 

52 

52 

52 

180° 

.605 

.999 

-2.2 

0.0 

.0555 

2 

47 

49 

5 

46 

48 

4 

49 

51 

50 

142° 

.549 

.909 

-2.61 

-0.42 

.0287 

5 

48 

52 

6 

45 

48 

7 

41 

40 

8 

55 

50 

9 

48 

47 

47 

116° 

.495 

.821 

-5.06 

-0.87 

.0276 

10 

48 

51 

11 

45 

47 

12 

45 

45 

» 

1? 

4l 

55 

14 

44 

47 

46 

101.50 

.475 

.784 

-5.26 

-1.06 

.0248 

15 

47 

49 

16 

44 

44 

17 

55 

55 

18 

47 

48 

19 

45 

47 

52 

78.5° 

.611 

1.010 

-2.15 

0.05 

.0519 

20 

46 

45 

21 

61 

62 

22 

5' 

56 

25 

5<* 

52 

24 

52 

52 

55 

64° 

.671 

1.122 

-1.74 

0.46 

.0576 

25  55  - 

26  61  59 

27  54  58 

28  54  58 


29 

— 

— 

60 

58° 

.797 

1.520 

-0.96 

1.21 

.0417 

50 

59 

64 

51 

57 

59 

52 

67 

67 

67 

0° 

1.000 

1.655 

0 

2.19 

.0556 

•Vertical  Average  -  The  averaged  output  of  all  beam*  of  equal  eleva- 
tion angle. 
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-6   t 


120° 


-4 


180°  '50* 

FIGURE  24    TRUE   FIELD  VERTICAL  DISTRIBUTION 
HERATIVE    SOLUTION 


52 


by  dashed  lines.  This  curve  is  used  to  select  values  of  X/     in  10 
degrees  steps  to  complete  the  numerical  integration  of  equation  (1) 
of  paragraph  6.1.  The  integration  of  aquation  (1)  was  accomplished 
with  incremental  steps  of  10  degrees  in  elevation  and  azimuth. 
Figure  25  represents  the  results  of  all  computations  with  the  ex- 
ception of  the  integration  over  the  variable,  9  .  The  ordinate  is 
G{<9)   X  (  8)  sin  9,   the  abscissa,  0   •  The  8  curves  represent  the  8 
vertical  beams  and  is  included  to  illustrate  the  fact  that  the  beam 
response  contributes  significantly  in  all  directions  and  emphasizes 
the  importance  to  this  process  of  knowing  the  beam  patterns  exactly. 
The  curves  of  Fig.  25  are  then  integrated  by  planimeter  and  the 
decibel  ratios  between  a  given  beam  and  the  beam  directed  at  the 
surface  obtained.  The  results  are  plotted  as  A"sin  Fig.  24.  The 
decibel  difference  between  dots  and  X's   is  measured  and  twice  that 
value  applied  as  a  correction  to  the  first  assumption,  dashed  curve, 
of  Fig.  24.  The  dotted  curve  results.  The  numerical  integration 
process  is  repeated  yielding  the  open  circles.  A  final  iteration 
yields  the  solid  curve  which  is  the  curve  of  Fig.  18. 


FIGURE  25.     G(0)X(0)SIN0  VS  8 


84 


BIBLIOGRAPHY 

1.  R.   A.   Fro«ch,    Unpublished  paper,   Acous.   Soc.   Am.  Meeting,   June 
(1960). 

2.  V.   0.   Knudsen,   R.   S.  Alford  and  J.   W.   Baling,   Rpt.  No.   S,   NDRC 
See.  No.  6.   1-NDRC-1848  (1944). 

5.  C.  Eckart,   J.   Acous.   Soc.   Am.   2£,    195  (1955). 

4.  R.  E.   Roberson,   J.   Acous.   Soc.  Am.   2£,   555  (1951). 

5.  R.  J.   Urick,   NRL  Report  5769   (1951). 

6.  C.  N.  Miller,   NEL  Conf.   Rpt.  698  of  22  June  1956. 

7.  NRDC  Summary  Technical   Reports,   Div.   6,   8. 

8.  C.  A.  Teer,   Great  Britain  UDE  Report  No.  65,    (1949). 

9.  V.  C.  Anderson,   J.  Acous.   Soc.  Am.  ^0,   440  (1958). 

10.  V.    C.   Anderson,   J.  Aoous.   Soc.   Am.  ^2,   867   (19^0). 

11.  P.   Rudnick,   J.   Aoous.   Soc.  Am.  ^2,   871    (i960). 

12.  H.  Primakoff,  Columbia  Div.   of  War  Res.  Library  of  Congress 
No.   L5#5 

15.     V.  C.   Anderson  MPL  Tech  Memo  112,   26  January  19^0. 

14.  R.  N.   Bracewell,   Aust.   Jour,   of  Phys.   7,   615  (1954). 

15.  T,  G.   Bell,    Paper  to  18th  Underwater  Acous.   Symp.,  Monterey, 
25  October  19^0. 

16.  K.   V.  MacKenaie,   J.   Acous.   Soc.   Am.  ^2,   221    (19&>). 

17.  H.  W.  Marsh  and  M.   Schulkin,   USL  Tech  Memo  No.   1110-110-54 
27  August  1954. 

18.  J.  L.  Lawson  and  G.   E.   Uhlenbeck,    "Threshold  Signals,"  McGraw- 
Hill,   M.I.T.   Radiation  Laboratory  Series,   Vol.   24. 


